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Abstract 

We calculate the two-loop renormalization group (RG) /3-function of a massless 
scalar field theory from the irreducible version of Polchinski's exact RG flow equa- 
tion. To obtain the correct two-loop result within this method, it is necessary to 
take the full momentum-dependence of the irreducible four-point vertex and the 
six-point vertex into account. Although the same calculation within the orthodox 
field theory method is less tedious, the flow equation method makes no assumptions 
about the renormalizability of the theory, and promises to be useful for perform- 
ing two-loop calculations for non-renormalizable condensed-matter systems. We pay 
particular attention to the problem of the field rescaling and the effect of the asso- 
ciated exponent r] on the RG flow. 
PACS numbers: 05.10.Cc, ll.lO.Gh, ll.lO.Hi 



1 Introduction 



There are two different ways of formulating renormalization group (RG) trans- 
formations: the first formulation [1], which we shall call the field theory 
method, relies on the fact that in a renormalizable field theory all physical 
quantities can be expressed in terms of a finite number of renormalized cou- 
plings which can be defined at an arbitrary scale /i. Because the bare theory 
is defined without reference to /i, any change in the renormalized correlation 
functions in response to a variation of /i must be compensated by a corre- 
sponding change in the renormalized couplings, such that the bare correlation 
functions remain independent of /i. This relation is usually expressed in terms 
of a partial differential equation for the renormalized correlation functions, 

^ Address after September 2000: Institut fiir Theoretische Physik, Johann Wolfgang 
Goethe-Universitat Frankfurt, Robert-Mayer-Strasse 10, 60054 Frankfurt am Main. 



Preprint submitted to Elsevier Preprint 



1 February 2008 



the Callan-Symanzik equation. The Gell-Mann-Low /^-functions describe how 
the renormahzcd couphngs must depend on ^ to guarantee that for fixed bare 
couphngs the bare theory is independent of /x. 

The other, more modern formulation of the RG was pioneered by Wilson 
[2]. In this approach, which has been very fruitful in statistical physics and 
in condensed matter theory [3], one starts from a local bare action S{(j)} 
with some ultraviolet cutoff Aq, and integrates out the degrees of freedom 
for momenta down to some smaller cutoff A < Aq. After this elimination 
of the degrees of freedom, the momenta and in general also the fields have 
to be properly rescaled in order to obtain a fixed point of the RG [4]. In 
this way one obtains an effective action with a new ultraviolet cutoff A. In 
general, the new action is a very complicated non-local functional of the new 
fields. The Wilsonian fiow equations describe how the various couplings of the 
effective action change as the scale A is reduced. An important advantage of 
the Wilsonian RG over the field theory method is that the Wilsonian method 
can also applied to non-renormalizable field theories. In fact, condensed matter 
systems are usually characterized by a finite lattice spacing, which plays the 
role of a physical ultraviolet cutoff. Hence, there is no need to perform the 
continuum limit in condensed matter theory. 

In the important work [5] Wegner and Houghton showed that the change 
of the effective action due to an infinitesimal change in A can be described 
in terms of a formally exact functional differential equation. However, for 
practical calculations the Wegner- Houghton equation has some unpleasant 
features, so that it has not been widely used to solve problems of physical 
interest. Recently, there has been a renewed interest in exact formulations 
of the Wilsonian RG. Polchinski [6] noticed that a particular form of the 
exact Wilsonian RG equation (which is now called the Polchinski equation, 
see Sec. 2.1) offers a straightforward approach to proof the renormalizabihty 
of quantum field theories. 

Due to its intuitive physical interpretation and its greater generality, the 
Wilsonian RG is the method of choice in condensed matter theory and statisti- 
cal physics [3] as long as a one-loop calculation is sufficient. In some problems, 
however, the correct behavior of physical observables requires the knowledge 
of the two-loop RG /5-function. An example is the asymptotic low-temperature 
behavior of the susceptibility of classical [7] and quantum [8] Hciscnbcrg mag- 
nets in two dimensions. However, till now even condensed matter physicists 
resort to the less intuitive field theory method for two- loop calculations [9] . In 
fact, the opinion seems to prevail that two-loop calculations within the exact 
Wilsonian RG are more difficult than within the field theory method. In this 
work wc shall show that this is not necessarily the case if one uses the irre- 
ducible version of the exact RG [10,11] in the form derived by Wetterich [12] 
and, independently, by Morris [13]. 
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We consider a simple scalar field theory with a bare action 



S{<P} = So{<P} + 5i„t{0} , 



(1.1) 



where the free part is 



(27r) 



0kG'o-^(k)0_k . 



(1.2) 



|k|<Ao 



Here (pk = J dre *'*''"0(r) are the Fourier components of a real scalar field 0(r) 
in D dimensions, and 



The interaction part Siat{4>} may also contain counter terms that are quadratic 
in the fields. Later wc shall assume that it is a local function of the fields that is 
invariant under (f) —0. The one-loop RG /3-function for this theory has been 
derived some time ago from the Polchinski equation [6] by Hughes and Liu [14]. 
Several strategies for calculating the RG /^-function at two loops have been 
proposed. A derivation of the two-loop /9-function entirely within the frame- 
work of the exact RG was given by Papenbrock and Wetterich [15], who used 
the exact evolution equation for the eff'ective average potential and worked 
with a smooth cutoff in momentum space. Because sharp cutoffs in momen- 
tum space generate long-range interactions in real space, Wilson and Kogut 
[2] recommended to use a smooth cutoff for calculations beyond one loop. On 
the other hand, the one-loop /^-function is most conveniently performed with 
the Wllsonian RG using a sharp cutoff, so that it would be useful to have a 
formulation of the Wilsonian RG where the same cutoff procedure can also 
be employed for higher order calculations. The precise relation between the 
orthodox field theory approach and the exact Wilsonian RG was established 
by Bonini et al. [16], who used this relation to obtain the /3-function at the 
two-loop order. Another strategy for calculating the two-loop /3-function of 
0^-theory was adopted Pernici and Raciti [17], who derived the Gell-Mann 
Low equation from the exact RG, and then used this equation to express the 
/5-function in terms of Wilsonian Green functions. In Ref. [13] an attempt has 
been made to calculate certain contributions to the two-loop coefficient of the 
/3-function by truncating the functional differential equations of the exact RG 
within a momentum scale expansion; unfortunately, this truncation is uncon- 
trolled (see also Ref. [18]) and the result is not a good approximation to the 
known two- loop result obtained within the field theory method [1]. Very re- 
cently Morris and Tighe [19] succeeded in calculating the two-loop /5-function 
from the exact RG by resumming the momentum scale expansion for sharp 
cutoff and the derivative expansion for smooth cutoff to infinite orders. 




(k) = k2 + m^ 



(1.3) 
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In this work we shall show how the two-loop /3-function can be obtained di- 
rectly within the framework of the exact Wilsonian RG using a sharp cutoff 
in momentum space, without resumming the momentum scale expansion, and 
without any reference to the Callan-Symanzik equation. Our motivation is 
to demonstrate that a two-loop calculation within the exact RG is concep- 
tually quite simple and involves only straightforward albeit tedious algebra. 
Our aim is to convince the reader that the exact RG is the most convenient 
method for performing two-loop RG calculations in condensed matter theory 
and statistical physics. 

The rest of this paper is organized as follows: We begin in Sec. 2 with a brief 
summary of the different versions of the exact Wilsonian RG. In Sec. 3 we 
discuss the proper rescaling of momenta and fields in the exact RG. In Sec. 4 
we explicitly write down the exact flow equations for the free energy, the irre- 
ducible two-point vertex, the irreducible four-point vertex, and the irreducible 
six-point vertex. We also elaborate on the relation between the anomalous 
dimension of the field and the irreducible two-point vertex. In Sec. 5 we make 
some general remarks concerning the structure of the RG flow. In Sec. 6 we 
briefly re-derive the one-loop /3-function. The two-loop calculation is then pre- 
sented in Sec. 7. Our conclusions can be found in Sec.8. Finally, in an Appendix 
we briefly summarize the definition of various types of generating functionals 
and introduce some convenient notations. 



2 Exact flow equations 



To integrate out degrees of freedom with momenta in the shell A < |k| < Aq, 
we introduce the cutoff-dependent free propagator with matrix elements 

in momentum space given by 

[<'^°]k,k' = {2nfS{k - k')<'^°(k) , (2.1) 



where 



Gr°(k) , zy-'--- . (2.2) 



^,(|k|,A)-^,(|k|,Ao) 



k2 + 



m, 



Here 9e{k, A) ^ 1 for A; > A, and 9^{k, A) ^ for /c < A, such that 9e{k, A) 
drops from unity to zero mainly in the interval |/c — A| < 2e. We are using 
here the same notation as Morris [13]. For explicit calculations, we shall later 
on focus on the sharp cutoff case, lim^^o ^) = ^(^ ~ but at this 
point, we shall keep e finite. There are several types of exact flow equations, 
corresponding to different types of correlation functions. A brief summary of 
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the various generating functionals and their representations in terms of func- 
tional integrals or functional differential operators is given in the Appendix. 
For convenience, let us now briefly summarize three types of exact flow equa- 
tions. Our discussion closely follows Refs. [20,21], but we shall be more careful 
to keep track of all field-independent contributions to obtain the fiow equa- 
tion for the free energy. Apart from the three flow equations given below, two 
more types of flow equations should be mentioned: The Wegner-Houghton 
equation [5], which was derived for system with discrete momenta and sharp 
cutoff, is closely related to the Polchinski equation [20]. Diagrammatically 
these equations contain terms which are not one-particle irreducible. If one 
works with a sharp cutoff in momentum space, this leads to ambiguities in the 
thermodynamic limit, where the momenta form a continuum. Nicoll, Chang 
and collaborators [10] avoided these difficulties by considering the flow equa- 
tion for the Legendre effective action, which generates irreducible correlation 
functions. The exact flow equation given in Sec. 2. 3 is equivalent to the flow 
equation of Ref. [11]. 



2. 1 Polchinski equation for Q^^° 



The generating functional of the amputated cutoff-connected Green functions 
^A,Ao pg^j^ deflned via (see Eq. (A. 10)) 



gSac'^° {</>} _ g|(^>G'o''^''^)g-S'int{</>} 



(2.3) 



Differentiating both sides with respect to A, we obtain the Polchinski equation 

[6], 



+ 



A,Ao "^ac 




5(t> 



, (2.4) 



where 5a on the right-hand side acts only on Gq'^°. Prom Eq. (2.3) we see 
that the initial condition at A = Aq is 



(2.5) 



Performing the functional derivatives in Eq. (2.4) in Fourier space, the Polchin- 
ski equation can also be written as 



dk 



^[aAGr°(k)] 



2 J (27r) 



(2.6) 
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This equation has been used in Ref. [14] to derive the one-loop /3-function for 
our simple scalar field theory. Let us emphasize again that diagrammatically 
the right-hand side of Eq. (2.6) contains terms that are one-particle reducible. 
This leads to technical difficulties in calculations beyond the leading order 
[13]. Although this problem can be avoided by expanding in terms 

of so-called Wick-ordcrcd monomials [23], for our two-loop calculation it will 
be more convenient to start from the exact flow equation for the generating 
functional of the irreducible correlation functions, see Eq. (2.12) below. 



2.2 Flow equation for Q^'^° 



Before discussing in some detail the (for our purpose) most convenient version 
of the exact RG, let us mention (for completeness) the flow equation for the 
generating functional Q^'^°{J} of the connected Green functions, which seems 
to be even less convenient than Eq. (2.4). From Eq. (A. 7) we see that ^^'^"{^} 
can be obtained from that of the corresponding generating functional of the 
amputated connected Green functions via 



(2.7) 



The exact flow equation is 



+ 



A,Ao 



-1^ 

SJ. 



tA,Ao 



-,9a(G'^^°) 



SJ 



-isg. 



A,Ao~ 



5J 



A,Ao 



(2.8) 



with initial condition at A = Aq 



^Ao,Ao|j| _ Q 



(2.9) 



This is an ill deflned initial problem, so that Eq. (2.8) is not very useful in 
practice. 



2.3 Flow equation for 



The exact RG equation for the the generating functional of the irreducible 
vertices turns out to be the most convenient starting point for our two-loop 
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calculation. Prom Eq. (A. 17) we have 



Gt'^^'W) = J) - 2 (Gr°)-V) - , (2.10) 



iA,Ao\-l, 



where J is defined as a functional of via 



5J{r) 



(2.11) 



The flow equation for Qtr^°{^p} is [12,13] 



A,AoN-l 



(2.12) 



where the functional V^'^°{</?} is the second functional derivative of ^i^'^°{</'}, 



(2.13) 



Unlike the authors of Refs. [12,13], we have kept track of the field-independent 
contributions to the right-hand side of Eq. (2.12), which determines the flow 
of the free energy [5]. Note that in the non-interacting limit, i.e. for V^'^° = 0, 
the right-hand side of Eq. (2.12) vanishes. The initial condition at A = Aq is 
simply 



(2.14) 



It is convenient to separate at this point from Eq. (2.13) the fleld-independent 
part, which can be identifled with the irreducible self-energy. 



(2.15) 



where by definition U^'^°{(p = 0} = 0. Then Eq. (2.12) can also be written as 



^^A,Ao^2^A,Ao ^2 _|_ (^A.,Ao^A,Ao~j ^ 



, (2.16) 
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where 



(2.17) 



is the exact two-point function. Note that in momentum space 
KPW} = (27r)^(5(k + k')E^'^o(k) , 



,, g,(|k|,A)-g,(|k|,Ao) 
k2 + mg + E^.^o(k) ■ 



(2.18) 
(2.19) 



For our two-loop calculation, it will be useful to take the sharp cutoff hmit of 
Eq. (2.16). In this limit 9e{k, A) 6{k — A). Performing the trace in momen- 
tum space, Eq. (2.16) reduces to [13] 



F) c^.^o__i [ ^(|k|-A) 

2] (27r)^A2 + mg + EAA,(k) 



-1 



k,-k 



^/^.m-A).n 



A2 + m2 + E^'^o(k) 



A2 + 



rrin 



(2.20) 



where V is the volume of the system. In deriving Eq. (2.20) we have ignored 
possible ambiguities that might arise from functions ^^(0) contained in the term 

Z^^'^o (1 + See Ref. [13] and the remark in Sec.4.3 for a further 

discussion of this point. 



3 Rescaling momenta and fields 



As it stands, Eq. (2.20) describes the change in the Legendre effective action 
due to the elimination of degrees of freedom with momenta in the interval 
[A, Aq]. A complete Wilsonian RG transformation includes also a rescaling of 
momenta and fields. A convenient way of including these rescaling transfor- 
mations into the above fiow equations it to rewrite them in terms of dimen- 
sionless variables by multiplying all quantities by appropriate powers of A. By 
taking the A-derivative while keeping the dimensionless variables constant, the 
rescalings are implicitly carried out [20]. In our case, it is useful to introduce 
dimensionless variables as follows, 

k = Aq , A = Aoe-* , (^k = A''''~''\/^<^q ■ (3.1) 

Here D'^ = ^[D — 2) is the scaling dimension of the field (p{r) in real space 
(so that D'^ — D is the scaling dimension of its Fourier transform (/7k). The 
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scale- dependent dimensionless factor Zt is the wave-function renormalization 
factor, which is related to the anomalous dimension 7]% of the field via 



r]t^-dMZt. (3.2) 



At this point the introduction of Z^, which cannot be derived from dimensional 
analysis, seems to be a rather ad hoc procedure. However, in general it is 
necessary to introduce such a factor, because otherwise it may happen that 
the RG transformation never reaches a fixed point [4,5,20]. Given Eq. (3.1), it 
is (almost) natural to define 

^ir'^n^k}-ylnZ, = ^,{(^J, (3.3) 

W^'^n^k} = A-^^"Zr^W,{^J , (3.4) 
Q^M ^ AS^^ZtGt , (3.5) 
+ + E^'^°(k) = A^Z^ndql) , (3.6) 

where have used the fact that E'^'^"(k) depends only on |k|. Here is the 
number of Fourier components in a system with ultraviolet cutoff A = AqC"*, 

Nt-y \ (^^(Aoe-* - |k|) = ^T/A^-^* ^ AToe-^* , (3.7) 
where Kj^ is the volume of the unit sphere in D dimensions divided by (27r)^, 

" 2^-%^/2r(L'/2) ■ ^^-^^ 

The constant ^^ In Zt in Eq. (3.3) modifies the scahng of the free energy. This 
term is due to the fact that the anomalous field rescahng leads also to a change 
in the integration measure, as explained in detail by Wegner and Houghton [5]. 
In our case, field rescaling 0k A^^~^Zt0k modifies the generating functional 
of the cutoff connected Green functions as 

^ , (3.9) 



so that 



gA,Ao ^gA,Ao + !^l^Zt, (3.10) 
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and hence, from Eq. (A. 17) 



(3.11) 



Substituting the above rescaled variables into Eq. (2.20) the left-hand-side 
becomes (after multiplication with a factor of A) 

AaAg^i^'^ny'k} = -dtGti^^} + 1ZsQt{^^} + ^[r]t + DlnZt] , (3.12) 



where TZs is the functional differential operator describing the rescaling of 
fields and momenta [5,20] 



(27r) 



r dc 
7(2^ 



~ 1 ^ 



(3.13) 



Here 



Dr^D^ + '^ = \{D-2-rm) 



(3.14) 



is the full dimension of the renormahzed field (in real space). Performing 
the angular integration in Eq. (2.20) in terms of D-dimensional spherical co- 
ordinates, we finally obtain the exact flow equation 



dtgt = nsQt + 



2rt(l) 



q,-q 



Nt 

+ -4 



r]t + Din 



(3.15) 



Here (. . .)q means averaging over all directions of the D-dimensional unit 
vector q. Note that Gt is a diagonal matrix in momentum space with matrix 
elements 



[G*U = (27r)^5(q + q')G,(q), 

'|q|-l)-^^(|q|-e*) 



(3.16) 
(3.17) 
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where rt{q) is defined in Eq. (3.6). The initial values ait — are by construc- 
tion 

G'i=o(q) = , Zt=o = l , n=o{q) = + (^)' ■ (3.18) 

Eq. (3.15) is the central result of this section. At a fixed point of the RG Qt 
become stationary, 

dtgt{<fci} = . (3.19) 

Of course, usually it is impossible to solve Eq. (3.15) exactly, so that we have 
to resort to approximations. 



4 Exact flow equations for the irreducible vertices 

To generate an expansion in the number of loops, let us we expand the func- 
tional Gti'fci} in powers of Fourier components (fq of the field, 

xri")(qi,...,q„)(^q,...^q„ . (4.1) 
The functional Ut on the right-hand side of Eq. (3.15) can be written as 

xri'^)(qi,...,q„)-i-— . (4.2) 

We now substitute Eqs.(4.1) and (4.2) into Eq. (3.15) and expand 

Ut(l + GtUt)-^ ^Ut- UtGtUt + UtGtUtGtUt - . . . . (4.3) 

Collecting all terms with the same powers of the fields on both sides of Eq. 
(3.15), we obtain an infinite hierarchy of flow equations for the irreducible n- 
point vertices r[^\qi, .... q„). If the bare action is invariant under (/? — > — (/?, 
only the vertices with even n are finite. It turns out that for the two-loop 
calculation of the RG /5-function it is sufficient to truncate the hierarchy by 
setting r^"'(qi, . . . , q„) = for n > 8. We now explicitly give the exact 
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flow equations for the irreducible vertices that are required for the two-loop 
calculation. 



4.1 Free energy 



The flow equation for vf'^ describes the flow of the interaction correction to 
the free energy. Prom the field-independent part of our exact fiow equation 
(3.15) we find 



rjt + Dln 



(4.4) 



The initial condition is 



^(0) 
t=0 



. 



(4.5) 



Setting r^''"' = Ntft and using Nt = Nqc'^^ (see Eq. (3.7)), wc obtain the 
flow equation for the interaction correction to the free energy per Fourier 
component, 



9t/t = i^/* + | + §ln 




(4.6) 



Note that in the corresponding flow equation for the potential Vq given by 
Wegner and Houghton (see Eq. (3.9) of Ref. [5]) the term rjt is replaced by 
T) — 2, where rj is assumed to be independent of t. The extra —2 is due to the 
fact that the potential vq introduced in Ref. [5] represents the total free energy, 
while our vf*^ contains only the interaction correction to the free energy. Hence 
for a non-interacting system our r|°^ does not flow and remains identically 
zero. The solution of Eq. (4.6) with the correct initial condition is 



rjt' + Din 



(4.7) 



4-2 Irreducible two-point vertex and self-energy 



From Eq. (3.15) we obtain the following exact flow equation for the dimen- 
sionless irreducible two-point vertex. 
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dtr?\ci, -q) = (2 - r?* - q • V^) rf ^(q, -q) 

/-p(4). 



The initial condition at i = is determined by the interaction part of the bare 
action Sint{<p}, see Eq. (2.14). If Sint{<p} does not have any term quadratic 
in the fields, then the initial condition is fI^q^c^, — q) = 0. However, to reach 
the Gaussian fixed point in D > 4 it is necessary to fine tune the bare action 
such that the system flows along the critical surface. As explained in detail in 
Sec.5, in this case the initial value t[%{0, 0) becomes a function of the bare 
four-point vertex (see Fig.l and Eq. (5.5)). Of course, such quadratic term 
in the fields can also be absorbed into a redefinition of the Gaussian part of 
the action, but we find it more convenient to include all interaction-dependent 
terms in 5'int{<^}- 

By definition, our dimensionless parameter rt{q) is related to the irreducible 
two-point vertex r[^\q,, — q) via 

n{q) = ZMQ) + ^?\^,-^) , (4.9) 



implying 



[a, - 2 + r;, + q • V^][r?\ci, -q) - n{q)] = 2Z,^ . (4.10) 



Hence, up to a re-definition of the four-point vertex by a momentum- 
independent constant proportional to ttIq/Aq, the fiow equations for rt{q) and 
rp-*(q, — q) are identical. Note that the authors of Ref. [11] prefer to work 
with rt{q) (their fiow equation for rt{q) is thus equivalent to our Eq. (4.8), 
except that these authors replace rjt by its fixed point value). In this work we 
find it more convenient to parameterize the two-point function in terms of its 
irreducible part r(^^(q, — q). 

Physically it is clear that the irreducible two-point function should uniquely 
determine the anomalous dimension rjt — —dt In Zf of the field. Hence, rjt on 
the right-hand side of Eq. (4.8) implicitly depends on r[^\q, — q). The precise 
relation becomes more transparent when we rewrite Eq. (4.8) in terms of the 
dimensionless irreducible self-energy at{q) defined by 

at{q) = A-^E^'^°(Aq) = Z^^rf (q, -q) , (4.11) 



where A = Aqc *. Using dtZ^ = —Ztrjt, it is easy to show that Eq. (4.8) is 
equivalent with 
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dtat{q)^{2-q-V^)atiq) 



+ 



2Z2[ro(l) + ai(l)] 



rf)(q,-q,q',-q': 



(4.12) 



The wave-function renormalization Zt on the right-hand side is now fixed by 
demanding that [22] 



(4.13) 



which imphes 



1 + 



dat{q) 



d{q^ 



(4.14) 



From the definition Eq. (4.11) it is then easy to show that Zt is related to the 
irreducible two-point vertex via 



1 - 



d{q^) 



(4.15) 



q2=0 



Using Eq. (4.14) to eliminate Z^ from Eq. (4.12), we obtain the following 
non- linear partial differential equation for (Jt(g), 



dtat{q) = {2-Ci-V^)at{q) 



K 



D 



2[ro(l) + a,(l)] 



d(yt{q) 



d{q^) 



rfHq,-q,q',-q')). -(4.16) 



Given a solution of Eq. (4.16) with appropriate initial condition, the anomalous 
dimension of the field can be obtained from Eq. (4.14). Note, however, that 
the condition (3.19) that the RG has a fixed point requires 



9trf)(q,-q) = 0, 



(4.17) 



which according to Eq. (4.11) means that 

dtCTtiq) ^Wt{q) , 



(4.18) 



where r] is the value of r]t at the RG fixed point. Thiis, for a system with a 
finite anomalous dimension the irreducible self-energy is not stationary at a 
fixed point of the RG. This is the reason why Eq. (4.8) is more useful than 
Eq. (4.12) if we are interested in possible fixed points of the RG. Keeping 
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in mind that dtO't{q) on the left-hand-side of Eq. (4.18) can be replaced by 
the right-hand side of Eq. (4.16), we see that Eq. (4.18) can also be viewed 
as a non-linear eigenvalue equation [20,24]. From this point of view it is not 
surprising that a RG fixed point can only be reached for certain discrete values 
of ?7 [4]. 



4-3 Irreducible four-point vertex 



The irreducible four-point vertex satisfies the following fiow equation. 



^trf^(qi,q2,q3,q4) 



4-D-277t-^q,-Vq, 

i=l 



rr^(qi,q2,q3,q4) 



K 



D 



2rt(l) 



rf^q^ -q, qi,q2,q3,q4) 



Kr 



D /■p(4) 



rt{l) 



n n4 -q-qi -q2,qi,q2)Gt(q + qi + q2) 



X rf ^(q + qi + q2, -q, qa, q4) + (q2 ^ qa) + (q2 ^ q4)). ■ (4.19) 

As pointed out by Morris [13], this equation is only valid if none of the com- 
binations of external momenta qi + q2, qi + qs, and qi -|- q4 that enter the 
propagators ^((q-l-qi -|-qj) vanishes. Otherwise these expressions contain the 
function ^(0), which is ambiguous. However, in physical problems where noth- 
ing drastic happens at vanishing momenta these special points in momentum 
space should not matter in the calculation of any observable. In this case we 
may simply ignore this problem, which we shall do from now on. On the other 
hand, if the field acquires a vacuum expectation value so that its q = Fourier 
component is finite, the sharp cutoff limit has to be taken more carefully, or 
the vacuum expectation value has to be treated separately [13]. 



4.4 Irreducible six-point vertex 



Eq. (3.15) implies the following flow equation for the irreducible six-point 
vertex. 



9*rf)(qi,...,q6) 



6-2L>- 377t-^q, • Vc 



rj'^(qi,---,q6) 



i=l 



K 



D 



15 



jy- 15 terms 

E (ri'Hq,-q-Qi,W(q + Qi)rr(q + Qi,-q,/2); 



nil) 



jy- 45 terms 

+^ E (rrnq,-q-Qi,W(q + Qi)rJ^)(q + Qi,-q + Q2,/3) 



{/l,/2},/3 



xGi(q-Q2)rJ'Hq-Q2,-q,/2)). ■ (4.20) 



We are using here the same notation as Morris [13]: = Sq^e/j qfc, and 
J2{ii,i2},ii,---im i^sans a sum over all disjoint subsets, Jj H /j = such that 
^/j = {qi, . . . q„}. The symbol {/i, I2} means that this pair is counted only 
once. The 15 terms of the sum involving the combination r'^'^^F*^^^ correspond 

to the 15 = possibilities of choosing distinct subsets of four momenta 

out of the six available momenta. The 45 terms in the last sum of Eq. (4.20) 

correspond to one half of the 90 ) possibilities of picking subsets of 



four momenta out of a set of six momenta, and then picking again subsets of 
two out of these four momenta. 



5 Relevant and irrelevant couplings 



Before embarking on the two-loop calculation of the RG /3-function, we should 
identify the relevant and irrelevant couplings. 



5.1 Relevant couplings 



In D > 4 the only fixed point of the RG flow is the Gaussian fixed point, 
where all irreducible n-point vertices with n > 2 vanish. Following the usual 
jargon, we call all couplings with a positive or vanishing scaling dimension 
relevant couplings. For D > 4 there are only two relevant couplings. The first 
is the momentum-independent part of the irreducible two-point vertex, 

//t = rf (0,0) = Ztat(0). (5.1) 



This coupling is strongly relevant and has scaling dimension +2. The second 
relevant coupling is the coefficient of the term of order in the expansion of 
Ti (q, — q) for small dimensionless wave- vectors q. 



Ct 



(q.-q) 
a(g2) 



(5.2) 
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see Eq. (4.15). The scaling dimension of q vanishes. Let us emphasize that in 
general q does not vanish at the Gaussian fixed point, but approaches a finite 
value which depends on the bare interaction, see Fig. 2 below. Instead of Ht and 
Ct, we could also parameterize the two relevant couplings in terms of (Tt(0) and 
Zt, but we prefer to work with nt and Cj, because these quantities are directly 
related to the derivatives of our dimcnsionless two-point vertex r(^'*(q, — q). 
In D < 4 the momentum-independent part of the four-point vertex 

gt = rl'\0, 0,0,0) (5.3) 



is also relevant. The scaling dimension of gt is positive in D < 4 and vanishes 
in D — A. All other couphngs are irrelevant in D > 4, so that the RG trajec- 
tory in the infinite-dimensional parameter space of all couplings rapidly fiows 
approaches the three-dimensional subspace spanned by the relevant couplings 
{Ht,Ct,gt}- As explained in detail by Polchinski [6], for sufficiently large t the 
irrelevant couplings can then be expanded in powers of the relevant couplings. 

From now on we shall consider the fiow on the critical surface in D > 4, which 
is the infinite-dimensional manifold in coupling space which fiows into the 
Gaussian fixed point under the RG transformation. Recall that in D < 4 the 
Gaussian fixed point becomes unstable, and a non-trivial fixed point emerges, 
the Wilson-Fisher fixed point [2,3]. For simplicity, we shall assume a vanishing 
bare mass (mo = 0), so that ro(g) = g^. The projection of the RG flow onto the 
plane spanned by gt and /it is shown in Fig.l. Obviously, at the Gaussian fixed 
point the mass renormalization /it vanishes, so that this fixed point represents 
a system at criticality. For sufficiently small gt the critical trajectory can be 
approximated by a straight line through the origin in the g-ii-plane, 

l^t = l^'^9t + 0{gl). (5.4) 

Note that Eq. (5.4) implies that, in order to reach the Gaussian fixed point, 
the initial value of the two-point vertex should be adjusted such that 

^xt=Q = rgo(o, 0) = 7^'^rgo(o, o, o, o) . (5.5) 

To determine the numerical constant 7*^^^ we use the fact that according to 
Eq. (4.8) 

dtpit^2fit + ^^^9t + 0{g^), (5.6) 

while Eq. (4.19) implies to leading order 

dtgt = {A-D)gt + 0{g^). (5.7) 
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K4gt 



Fig. 1. One-loop RG-flow for Z? = 4 in the (7-^-planc. Defining gt = Kj^gt, 
the one-loop flow is determined by dtfit = "^l^t + \gt/0- + f^t) and 
dtgt = {-i — D)gt — ^g]: / {1 + iit)"^ . The thick black line is the linear approxima- 
tion for the critical trajectory (/xt = —\gti see Eq. (5.4)), which flows towards the 
Gaussian fixed point (black dot). To reach this fixed point, the initial values gt=Q 
and iit=o must be fine tuned to lie on the critical trajectory. 

Here we have anticipated that the lowest corrections to r]t and Ct — 1 — Zt 
involve at least two powers of gt. Keeping in mind that ro(l) = 1, substituting 
Eq. (5.4) into Eq. (5.6) and using Eq. (5.7), we find 

' 2(D-2) ^ ' 



For the two-loop calculation of the /3-function presented below the linear ac- 
curacy given in Eq. (5.4) turns out to be sufficient. 

Having adjusted the strongly relevant coupling such that we reach the 
Gaussian fixed point, we are left with two relevant couplings, gt and q. Of 
course, in addition we have infinitely many irrelevant couplings, which are 
given by the higher order momentum-dependence of the two-point and four- 
point vertices, and the higher irreducible vertices. However, the irrelevant 
couplings can be expanded in powers of the relevant couplings, which is the 
key to perform a two-loop calculation within the fiow equation formalism. 
Hence, in general we expect that the projection of the exact RG fiow onto the 
subspace spanned by couplings gt and q can be described by equations of the 
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form 



dtCt = A{gt, ct) , dtQt = B{gt, q) , 



(5.9) 



with some dimensionless functions A{g,c) and B{g,c). 
5.2 Irrelevant couplings 

Because the RG trajectory rapidly approaches the manifold spanned by the 
relevant couphngs, all couplings with negative scaling dimension (the irrele- 
vant couphngs) become local functions of the relevant couphngs, which are 
independent of the initial conditions at t = 0. This is explained in detail in 
the seminal paper by Polchinski [6]. From the exact flow equation (3.15) for 
the irreducible vertices it is easy to see that for small gt and large t the leading 
behavior of the irrelevant coupling functions is 



Here the 7^"-'(qi, • • • , q™) are t-independent functions. Note that the leading 
corrections are independent of q. This is due to the fact that the leading term 
on the right-hand side of the flow equation (5.9) for Ct turns out to be propor- 
tional to g^, see Eq. (7.25) below. For a two-loop calculation within the exact 
RG we need to know the functions 7^^^(qi, q2, qs, q4) and 7^^^(qi, . . . , qe)- 



6 One-loop /9-function 

The one-loop /3-function has been calculated from Polchinski's flow equation 
by Hughes and Liu [14]. To this order, we may truncate the system of exact 

(n) 

flow equation by setting Tl = for n > 6. Furthermore, for a one- loop 
calculation the momentum-dependence of all couplings can be ignored. In 
particular, the wave-function renormalization can be neglected at one-loop 
order, so that Zt — 1 and hence Ct — rjt — 0. Obviously, the function A{gt,Ct) 
in Eq. (5.9) vanishes to this order. From the exact flow equation (4.19) we 
obtain for the momentum-independent part of the four-point vertex. 



rP(q, -q) - - c,q' ~ gh^'\q) + 0{g',) . 

rr^(qi, q2, qs, q4) -gt^ gh^^H^ii, q2, qs, q4) + o{g^) 

rJ^Hqi, ■■■:<ln)- 9:^W''\<lU • • • , qn) + 0(#+^)/^) . 



(5.11) 



(5.10) 



(5.12) 



dtgt = (4 - D)g^ - 3Kn (Gt(q))q^*' + 0{gl c^gl) , 



(6.1) 
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where Gt(q) is defined in Eq. (3.17). Note that by definition Gt=o{<i) = 0) but 
for t > and |q| < 1 the uhraviolet cutoff ^^(|q| — e*) in Eq. (3.17) can be 
neglected, so that to leading order we may approximate 



Using 

'1 2 ' 



(C(q)), = I , (6.3) 



we finally obtain for the right-hand sides A{gt,Ct) and B{gt,Ct) of the fiow 
equations (5.9) for the relevant couplings 

A{g„ct)^0{g^), (6.4) 
B{gt, ct) = (4 - D)g, - ^-Kj,gl + 0{gl) . (6.5) 

Eq. (6.5) is the known one-loop /3-function of massless 0^-theory [1]. Recall 
that in four dimensions = 2/(47r)^. 



7 Two-loop /3-function 

To obtain the flow equation for g^ at the two-loop order, we substitute 
Eqs.(5.11) and (5.12) into the exact flow equation (4.19) for the irreducible 
four-point vertex and collect all terms up to order gf. After a straightforward 
calculation we obtain for the generalized /9-function defined in Eq. (5.9), 

B{gt, Ci) = (4 - L> - 2r)t)gt ~ ^Kog^ 

+^Kngl [l^'^ - (T^'Hq, -q, 0, 0))^ + ^ {^^''\ci, -q, 0, 0, 0, 0)). 
+0{gt) . ■ (7.1) 

Recall that by definition 

77* = -a,ln(l-Q) = T^, (7.2) 

J- — Ct 

so that to this order the function B{gt,Ct) depends also on the second rele- 
vant couphng Ct- The term 7^^) in the square braces of Eq. (7.1) is due to the 
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interaction-dependence of the mass renormalization fit on the critical trajec- 
tory, see the thick sohd arrow in Fig.l and Eq. (5.4). The numerical value of 
^(1) is given in Eq. (5.8). The second term in the square brace of Eq. (7.1) 
is due to the momentum-dependent part of the four-point vertex, while the 
last term is due to the six-point vertex. Hence, to calculate the RG flow at 
the two-loop order, we have to calculate the anomalous dimension to order 
gf, the momentum-dependent part of the four-point vertex to order gf, and 
the six-point vertex to order gf. Note that it would be incorrect to ignore the 
momentum-dependence of the six-point vertex, because in Eq. (7.1) we need 
rf)(q, -q, 0,0, 0,0) for |q| = 1. 

7.1 Momentum-dependent part of four-point vertex 

Let us begin with the calculation of the momentum-dependent part of the 
four-point vertex, 

ff ^(qi, q2, qs, q4) = ^t^\(li, q2, qa, q4) - 9t ■ (7.3) 

Assuming that the bare four-point vertex is momentum-independent, 
(qi, q2, qa, q4) satisfies the initial condition at t — 0, 

f2o(qi,q2,q3,q4) = o. (7.4) 

Subtracting from the exact flow equation (4.19) for the full four-point vertex 
the corresponding equation with all external momenta set equal to zero, it is 
easy to show that to order g^ the function f ^'^-'(qi, q2, qa, q4) satisfies 

ai + D-4 + ^q,-Vq, 

1=1 

where 

q2, q3, q4) = -Kd{C{c^ + qi + qa) + C(q' + qi + qa) 

+C(q' + qi + q4)-3C(q'))., . (7.6) 

On the right-hand side of Eq. (7.6) we have replaced the exact propagators 
G't(q) by the non-interacting propagators without ultraviolet cutoff C(q), see 
Eq. (6.2). Such a replacement is only justified for t ^ 1, so that strictly speak- 
ing the approximation (7.6) cannot be used to study the fiow of Eq. (7.5) 
in the interval < t ^1. However, as discussed by Polchinski [6], the flow 



rj^^qi^ q2, qs, q4) = /^'^Hqi> q2, qa, q4) , (7.5) 
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of irrelevant couplings for large t becomes independent of the precise initial 
conditions (see below), so that for our purpose of investigating the infrared 
properties of the system we may ignore this subtlety. 

To solve Eq. (7.5) with initial condition (7.4), it is convenient to introduce the 
auxiliary functions 

f 2(qi, q2, qs, q4) = rl^^(e--^qi, e-^qs, e-^qs, e"^q4) , (7.7) 

^'Hqi, q2, qs, q4) = /^'He-^qi, e'^qs, e-^qs, e'^c^,) . (7.8) 

Then the solution of Eq. (7.5) is given by 

ff ^(qi, q2, qs, q4) = fJJ^o(qi, qa, qs, q4) , (7.9) 

where 

[d^-ds + D- 4] fJ5(qi, q2, qs, q4) = 9llf\^i. q2, qs, q4) • (7.10) 

The initial condition (7.4) implies 

f£o,.(qi,q2,qs,q4) = 0. (7.11) 

Because in Eq. (7.10) the external momenta q^ appear simply as parameters, 
this equation can be considered as a first order partial differential equation 
with two variables, t and s. The solution of Eq. (7.10) with the initial condition 
(7.11) is easily obtained, 

t 

fJ5(qi, q2, qs, q4) = / drgle-^^-'^^'-^^ q2, qs, q4) ■ (7.12) 



Setting s = in Eq. (7.12), we find after a shift of the integration variable 

t 

f J'^ (qi, q2, qs, q4) = / drgl^e-^^-'^^ if^ (q^, q^, qs, q4) ■ (7.13) 



At the first sight it seems that vf'\(^i, q2, qs, q4) is a non-local function of gt. 
However, by construction the term e^^^^'^'>'^ lf\c\x^ q2, qs, q4) in the integrand 
vanishes for large r at least as fast as e~'^^~^^'^ ^ so that the r- integration is 
effectively cut offatT^l/(L>— 3). For i^l/(D— 3)we may then approximate 
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g^_^ ^ in the integrand of Eq. (7.13), so that we finally arrive at the behavior 
anticipated in Eq. (5.11), 

fS'Hqi,q2,q3,q4)~t?tV'Hqi,q2,q3,q4) , t>l/{D-3), (7.14) 

with 

oo 

q2, qs, q4) = / o?Te-(^-^)V(^)(qi, q^, qa, q4) • (7.15) 



According to Eq. (7.1), for the two-loop flow equation for we need 

1 

^(4) ^ (^(4)(4_q,o,0)). = / (iAA^-^(/W(Aq,-Aq,0,0))^ , (7.16) 



where we have substituted A = e~'^. Note that the value of 7^^^ is determined 
by all momentum scales up to the infrared cutoff A. The momentum scale 
expansion discussed by Morris [13,18] corresponds to expanding the integrand 
in Eq. (7.16) to some finite order in A. Because the upper limit for the A- 
integration is A = 1, this expansion is not controlled by a small parameter 
[13,18]. Obviously, the correct numerical value of the two- loop coefficient of the 
/3-function can only be obtained if all terms in the momentum scale expansion 
are resummed [19]. Using the expression for /'^'^^(qi, q2, qs, q4) given in Eq. 
(7.6) and introducing D-dimensional spherical coordinates with cost? = q - q', 
we find 



iljy J J 



X 



^(A + 2cost9) 
l + 2Acos'i? + A2 



-^(cost?) 



(7.17) 



where — {2t:)^Kd is the volume of the unit sphere in D dimensions. 
Rearranging terms, 7*^^^ can also be written as 



V2 1 

j d-dism^f-^ J dXj 



A4-^(2cost? + A) 
T2Xcos¥+JF 



7r/6 1 

- J d^icos^f-^ J dA- 



2sini? 



2Asini? + A2 



(7.18) 
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where in the second integral we have shifted '»?—>■'»? — 7r/2. In Sec.7.5 we shall 
explicitly evaluate this integral for D — i. 



7.2 Anomalous dimension 

Having determined the momentum-dependent part of the four-point vertex 
to order g^, we can calculate the momentum-dependent part of the two-point 
vertex to the same order. In analogy with Eq. (7.3), let us define 

ff)(g)^rf)(q,-q)-/.,. (7.19) 

Note that according to Eq. (5.10) we expect for large t, 

r?\q) - + gH'\q) + 0{g^) . (7.20) 

Hence, the flow of our second relevant coupling q (which according to Eq. (7.2) 
determines the anomalous dimension rjt) can be obtained from the quadratic 
term in the expansion of Ti^\q) for small q. The function 7^^^(g), which by 
construction vanishes faster than q"^ for small q, is not needed for a two-loop 
calculation. 

Subtracting from the exact flow equation (4.8) for the irreducible two-point 
vertex the same equation with q set equal to zero, we obtain to leading order 
in gt 

[dt-2 + q. r?\q) = g^I^'\q) + 0{g^) , (7.21) 

where 

i^'Hq) = ^ (T^'Hq, -q, q', -q') - t^'Ho, o, q', -q'))., ■ (7.22) 

Expanding both sides of Eq. (7.21) in powers of q, we obtain to leading order 

dtct = a2g^ + 0{g^) , (7.23) 



where 



q;2 



9/(2) (g) 



d{q') 



(7.24) 



^2=0 
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From Eq. (7.23) wc conclude that the weak couphng expansion of the function 
A{gt, Cf) defined in Eq. (5.9) is 

A{gt, ct) = a2g^ + 0{gl gtCt) . (7.25) 
Eq. (7.23) is easily integrated, 

t 

Ct = a2 j dt'gl . (7.26) 



Hence, the relation between q and gt is non-local, so that in general also the 
anomalous dimension rjt given in Eq. (7.2) is a non-local function of gt. This 
is due to the fact that q is not irrelevant. However, to leading order the term 
in the denominator of Eq. (7.2) can be neglected, so that 

m = dtct + 0{gt) = a^gl + 0{gl) , (7.27) 



i.e. to this order rjt is a local function of gt- 

For an explicit calculation of the number 0:2, we need to calculate the leading 
gr-dependence of the function I^'^\q). Using Eqs.(7.22), (7.15) and (7.6), we 
find 

1 

= -Kl j dXX'-^ (/(A|q' + q|) - /(A))-, , (7.28) 




where we have defined 



.,^._/ ^(|q^^ + Ae|)-l \ _nn-i f {I - x^)^ 

^^^^ - \ (q" + Ae)2 /.„ - ^^l + 2Ax + P • ^^-^^^ 



Here e is an arbitrary constant unit vector. Apart from a different nor- 
malization, our function /(A) agrees with the corresponding function intro- 
duced in Appendix D of Ref. [11], were the anomalous dimension at the 
Wilson-Fisher fixed point has been calculated to order (4 — D)'^. Expanding 
/(2)(g) = a2q^ + 0{q^) we find from Eq. (7.28) 



a2 



2D 



1 _L 

(D-l) J dAA^-^/'(A) + J (iAA^-^/"(A) 



(7.30) 



In Sec.7.5 we shall explicitly evaluate Eq. (7.30) in four dimensions. 
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7.3 six-point vertex 

To complete the two- loop calculation of the function B{gt, ct) given in Eq. 
(7.1), we need to know the number 

7^'^ = (7^'nq,-q,o,o,o,o)). . (7.31) 

Recall that according to Eq. (5.12) the dimensionless function 7^^^(qi, . . . , qe) 
is defined in terms of the large f-behavior of the irreducible six-point vertex, 

(qi, • • • , qn) 5i*7^'Hqi, • • • , q^) + 0{gt) . (7.32) 

From the exact flow-equation (4.20) we find that to order gl the irreducible 
six-point vertex satisfies 

6 

at + 2D-6 + ^qi- Vq, 

i=l 

where 

45 terms 

7(«)(qi,...,q6)=X,, ^ (C(q' + Qi)C(q'-Q2))., . (7.34) 

{h,l2},h 

Eq. (7.33) can be solved with the same method as Eq. (7.5). The solution with 
initial condition r^f,^o(qi) ■ ■ ■ > qe) = is 

rf ^(qi, . . . , qe) = / dTgUe<'^-'^n^'\e-^c^,, . . . , e^q^) , (7.35) 



which is the analog of Eq. (7.13). For t ^ 1/(21^ — 6) wc may approximate 
g^_^ Ri g^ under the integral, so that Eq. (7.35) indeed reduces to Eq. (5.12), 
with 

1 

7^'Hqi, • • • , qe) = / rfAA^^-^/(^)(Aqi, . . . , Aqg) , (7.36) 



where we have substituted A = e~'^ . For our two- loop calculation we only need 
7(6) (q, -q, 0, 0, 0, 0). Using Eq. (7.34) this can be written as 



rf (qi,...,q6) = ^i*/^'nqi,---,q6), (7.33) 
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1 

7('nq,-q,o,o,o,o) = i^^ J dw^""-' 



+ 12 (C(q' - Aq)C(q') + C(q' + Aq)C(q') + + Aq))., 



(7.37) 



For the evaluation of (C'^(q'))^, one should be careful to properly define the 
value of the step function 6{x) at x = 0. Following the procedure discussed by 
Morris [13], we find 



(7.38) 



To obtain the number 7''^'* defined in Eq. (7.31), wc should average Eq. (7.37) 
over all directions of the unit vector q. Using D-dimensional spherical coordi- 
nates we obtain 



7 



(6) 



2{D - 3) 
xe(A + 2cosT?) 



+ 



1 



l + 2Acos'(? + A2 (l + 2Acos'(? + A2)2 
For an exphcit evaluation of Eq. (7.39) in D = 4 see Sec.7.5. 



(7.39) 



7.4 P -function 



Collecting all contributions to the function B{gt, q) given in Eq. (7.1) and us- 
ing the fact that the dependence on q enters only via the anomalous dimension 
Vt — C(2gt + 0{gl), we see that at the two-loop order the flow equation (5.9) 
for Qt can be written as 



(7.40) 



with 



/3o = 4 - L> , /3i = 



3K 



D 



, P2 = + PI + /3f + /3f , (7.41) 



where 



2(D - 2) ' 



(7.42) 
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/32' = -2a2, (7.43) 
(3i'^^-3Ku^^'\ (7.44) 

/3f = ^7^r (7.45) 

The term arises from flow of the momentum-independent part of the two- 
point vertex along the critical surface, see Eq. (5.8). The second term, (32, is 
due to the flow of the anomalous dimension, /^g^"* is due to the momentum- 
dependent part of the four-point vertex, and (32 is due to the six-point vertex. 
Note that the numerical values of (3!^^ and (32^^ are determined by the behavior 
of the four-point and six-point vertices at finite momenta, and thus contain 
the effect of an infinite number of irrelevant couplings. 



7.5 Four dimensions 



We shall now explicitly evaluate the two-loop coefficient (32 in D — 4, and 
show that it agrees with the result obtained within the field theory approach 

[1]- 

According to Eq. (7.42), the term (32 is in D = 4 given by 

Next, consider the contribution 02 due to the anomalous dimension. After an 
integration by parts we obtain from Eq. (7.30) in D = 4, 

«2 = ^ [2/(1) - 2/(0) + /'(I)] . (7.47) 
The required numbers are (see also Appendix D of Ref.[ll]) 



so that 



«2 = ^ . (7.49) 



Thus, the fiowing anomalous dimension in D = 4 is 



Vt = + 0{gl) , (7.50) 
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in agreement with the corresponding result of the field theory method [1]. Prom 
Eq. (7.43) we conclude that the contribution to the anomalous dimension of 
the field to the two-loop coefficient of the /^-function is 

= ■ (7.51) 



The calculation of P2 — — 3X^7^^^ is quite tedious. We first perform the 
A-integrations in the expression for 7*^^^ given in Eq. (7.18). In D — A the 
integrals can be evaluated analytically, 





1 

/ 



dX 



2sini? 



A(l-2Asin^? + A2) 



TT 
4 



1 



tani? ln[8(l — sini?)] — In sin 



(7.53) 



Substituting Eqs.(7.52) and (7.53) into Eq. (7.18), most of the t^-integrations 
can also be performed analytically. After several integrations by parts the final 
result for (3^^^ can be cast into the following form 



(4) 



K 



1 3\/3 ^ 6 /TT 
- + — 21n2 + -L - 

2 27r TT V3 



(7.54) 



where 



L{x) 



j (i'i? In cost? 




(7.55) 



is known as Lobachevskiy's function [25]. With L{^) 
merical value of jS'^'' turns out to be 



0.218391 ... the nu- 



(4) 



-Ki X 0.642204 . . . 



(7.56) 



Finally, let us evaluate the contribution — ^^^^^ from the six-point 
vertex to the two-loop coefficient of the /3-function in D = 4. To calculate the 
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integral 7^^^ in Eq. (7.39), it is convenient to first integrate by parts. Writing 

1 1 d 1 



(l + 2Acos'(? + A2)2 2Asini?ai?l + 2Acos'(? + A2 ' 
we obtain after some rearrangements 



^(6) ,^2(7 3^3 6 /• ^ 2Asin2^9-cos7? 







2Acosi? + A2 



(7.57) 



t/6 1 „ 

6 /■ ,0 f 2Acos^i? + sinT? I 

H — d-d / dX ^ — } . (7.58) 

TT J J l-2Asini? + A2| ^ ' 



Tlie A-integrations can now be performed analytically, 



/ d\ ^ 7^7 = coti? + sin 2^9 +sin2't?ln 2 1 + cos^? , 

J 1 + 2Acosi? + A2 2 ^ ^ ^ ^ 







l + 2Acos'(? + A2 

(7.59) 



1 

/ 



dA 



2 sin I 



2A cos^ -i? + sin-i? 
l-2Asin^? + A2 

fn 3'd 



[tani? + sin(2i?)] + cos^i?ln[2(l — sin^ 



(7.60) 



Substituting Eqs.(7.59) and (7.60) into Eq. (7.58) and performing the 
integrations we obtain after some tedious re-shufflings 



11 3v^ ^ 6/71 
— + 21n2 L - 

4 27r TT V3 



(7.61) 



In deriving this expression, we have used the fact that for < a; < ^ 
Lobachevskiy's function satisfies the functional relationship [25] 



L{x) 



TT 



In 2 



2 V2 



2x 



(7.62) 



which implies in particular 



1 n ^ 6 , /7r\ 9 , fir 
-ln2 L - + -L - 

2 TT V3y TT V6 



. 



(7.63) 
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We conclude that the contribution of the six-point vertex to the two-loop 
coefficient of the /^-function has the numerical value 

/3f ) ^Klx 2.892204 .... (7.64) 



Neither (3^ '/Kl nor ^r/Kj is a rational number, but their sum is: adding 
Eqs.(7.54) and (7.61), we find 

Pi'^+Pi'^ = lK!. (7.65) 



Collecting all terms, we obtain for the two-loop coefficient of the RG /3-function 



?(6) 



3 



1 

12 



Kl = ^Kl , (7.66) 



which agrees with the result obtained by means of the field theory approach 
[1]. The two-loop RG-flow in the plane of the two relevant couplings gt and q 
in D — 4: is shown in Fig. 2. 



8 Conclusions 



In this work we have shown how the known two-loop results for the RG (3- 
function of a simple massless scalar field theory can be obtained from the 
exact Wilsonian RG. Although from a technical point of view the two-loop 
calculation within the exact RG is more tedious than the same calculation 
within the orthodox field theory method [1], the Wilsonian RG is conceptually 
simpler, because there is no need to invoke rather abstract concepts such 
as dimensional rcgularization or minimal subtraction. Moreover, the exact 
RG has the advantage that it can be used even if the problem of interest 
cannot be mapped onto a renormalizable field theory. This is of particular 
importance in condensed matter systems, where an underlying lattice spacing 
always furnishes a physical ultraviolet cutoff, and there is no need for taking 
the continuum limit. An interesting system where the orthodox field theory 
approach fails is a two-dimensional electronic system with a Fermi surface 
that contains saddle points [26]. These give rise to van Hove singularities in 
the electronic density of states, which in turn generate non-renormalizable 
singularities in perturbation theory. As a consequence, this problem cannot 
be mapped onto a renormalizable field theory [27] . 

Due to its generality and conceptual simplicity, the exact Wilsonian RG is 
very popular in condensed matter theory and statistical physics. However, 



31 







0.5 1 

K^gt 



1.5 



Fig. 2. Two-loop RG-flow of the masslcss theory for D = i'm the g-c plane. The flow 
is determined by dtCt = jggt and dtgt = —^gt + li^t > where gt = K^gt. The thick 
line represents the Gaussian fixed point manifold, which can be parameterized by 
the value of the wave-function renormalization Zt=oo = 1 ~ Q=oo5 or alternatively, 
the initial interaction gt=o- Note that in Z? = 4 the two- loop /3- function (7.40) 
vanishes at = ^ « 1.058823 .... However, this is not a fixed point of the RG, 

because for gt > g* the second relevant coupling q flows to strong coupling. Because 
this graph is based on an expansion of the RG flow equations in powers of gt and 
ct, only the weak coupling regime K^gt <C 1 and q <C 1 can be trusted. 



even condensed matter theorists use the field theory method for performing 
two- loop RG calculations [9] . In the last few years the exact RG has been used 
to study the fermionic many-body problem [28] . To the best of our knowledge, 
two-loop RG calculations have not been performed for fermionic systems in 
D > 1. But also for one-dimensional fermionic systems, where the conven- 
tional field theory version of the RG method has been applied for many years 
[29], an ambiguity inherent in the conventional way of performing two- loop 
calculations has recently been pointed out [30]. It seems to us that the exact 
Wilsonian RG, in its irreducible form used in this work, is the most promising 
method for performing two-loop calculations for fermionic condensed matter 
systems. This formulation of the exact RG can be used even if the continuum 
limit does not exist, and has the potential of clarifying ambiguities of the type 
discussed in Ref. [30]. 
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Appendix: Four types of generating functionals 

Here we briefly summarize some basic definitions and representations of gen- 
erating functionals for various types of correlation functions. Although most 
of the following equations can be found in textbooks [1], we believe that this 
Appendix greatly enhances the readability of this work by giving a concise 
summary of the relevant expressions. 



A.l Disconnected Green functions 

The generating functional of the disconnected Green functions is given by the 
normalized partition function in the presence of external sources, 



G{J} 



Z{J} 



(A.l) 



where we have used the notation 




(A.2) 



The disconnected n-point functions are defined by 



r„) = (0(ri) ...(/)(r„)) 



(A.3) 



5J(ri)...5J(r„) 
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A. 2 Connected Green functions 



The generating functional Gc{J} of the connected Green functions can be 
defined by 



Note that by definition —Qc{J = 0} = F—Fq, where F is the free energy of the 
interacting system, and Fq = — In Zq is the free energy of the non- interacting 
one. The connected n-point functions can then be written as 



G(")(ri,...,r„) = ((/)(rO...0(r„)), 



SJ{r,)...SJ{rn) 



(A.5) 



j=o 



In contrast to Eq. (A.l), we have normaUzcd the integral in Eq. (A. 4) with the 
non-interacting partition function, such that — ^c{0} is just the change of the 
free energy due to the interaction. This normalization is useful, because then 
Qc{J} can be written in terms of functional differential operators as follows. 



= 62 



-Sir^t{4>} 



=GoJ 



(A.6) 



A. 3 Amputated connected Green functions 



The generating functional of the amputated connected Green functions is 



Gac{4>} 



Go{J}-^{J,GoJ) 



^c{Go-V}-^(0,Go^ 



(A.7) 



The amputated connected Green functions are then 



Gi:)(ri,...,r„) = 



5(/)(ri) . . . 50(r„) 



(A.8) 
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From the definition (A. 7) and Eqs.(A.4) and (A.l), it is easy to show that 
Qac{4>} has the following functional integral representation, 

JV{(f)'}e-s<^i^'} ' ^ ^ 



Furthermore, from the definition (A. 7) and the second hne of Eq. (A. 6) we see 
that Qa,c{(f)} has the following representation in terms of a functional differen- 
tial operator, 

g6ac{<?i>} _ Q^^M'^Oj^) ^-Sint{<t>} _ (A. 10) 



A. 4 Irreducible Green function 



The generating functional of the one-particle irreducible Green functions con- 
tains the information about the many-body correlations in the most compact 
form. This generating functional is obtained from the generating functional of 
the connected Green functions via a Legendre transformation. 



(A.11) 



where J{(fi} is defined as a functional of the classical field (f via 



6J{r) 



(A.12) 



Hence, 



00 1 

£M = F-Fo+5:- / dr,... dr„L(")(ri,...,r„)^(ri)...^(r„) , 

n=l •' 

(A.13) 



where 



LH(ri,...,r„) 



(A.14) 



Note that in the absence of interactions {Si^t = 0) the only non-zero vertex is 



L?Vi, r2) = Go ^(ri - = S{r, - r^) l-Vl + ml 



(A.15) 



35 



For the interacting system we have by definition 

L(2)(ri, rs) = G,\r, - r^) + S(ri - V2) , (A.16) 

where S(r) is the irreducible self-energy. To obtain the generating functional of 
all irreducible ra-point functions (including the two-point function), we should 
therefore subtract the free propagator from Eq. (A. 11), defining 



= C{^} - Go V) = iv, J) - liv, Gq-V) - Gc{J{v}} 



The irreducible n-point functions are then 



(A.17) 



Cir ^(ri, • • • , fn) 



(A.18) 



1/3=0 
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